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ABSTRACT 

The radiation by relativistic plasma particles is beamed in the direction of 
field line tangents in the corotating frame, but in an inertial frame it is aberrated 
toward the direction of rotation. We have revised the relation of aberration 
phase shift by taking into account of the colatitude of emission spot and the 
plasma rotation velocity. In the limit of small angle approximation, aberration 
phase shift becomes independent of the inclination angle a and the sight line 
impact angle (3. However, at larger altitudes or larger rotation phases, the shift 
does depend on a and f3. We have given an expression for the phase shift in the 
intensity profile by taking into account of aberration, retardation and polar cap 
currents. We have re-estimated the emission heights of the six classical pulsars, 
and analyzed the profile of a millisecond pulsar PSR J0437-4715 at 1440 MHz by 
fitting the Gaussians to pulse components. By this procedure we have been able 
to identify 11 emission components of PSR J0437-4715. We propose that they 
form a emission beam with 5 nested cones centered on the core. Using the phase 
location of component peaks, we have estimated the relativistic phase shift and 
the emission height of conal components. We find some of the components are 
emitted from the altitudes as high as 23 percent of light cylinder radius. 

Subject headings: pulsars: general — radiation mechanisms: nonthermal — stars: 
magnetic fields — pulsars: individual (PSR J0437-4715) 



1. INTRODUCTION 

The profile morphology and polarization of pulsars has been widely attempted to inter- 
pret in terms of emission in dipolar magnetic field lines (e.g., Radhakrishan & Cooke 1969; 
Sturrock 1971; Ruderman & Sutherland 1975; Lyne &; Manchester 1988; Blaskiewicz et al. 
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1991; Rankin 1983a&b, 1990, 1993; Hibschman & Arons 2001). Most of the radio emission 
models assume (1) radiation is emitted by the relativistic secondary pair plasma, (2) beamed 
radio waves are emitted in the direction of field line tangents, and (3) emitted radiation is 
polarized in the plane of dipolar field lines or in the perpendicular directions. 

Prom the theoretical point of view, it is highly preferable to know the precise altitude 
of radio emission region in the pulsar magnctosphere. By knowing the emission altitude, 
one can infer the probable plasma density, rotation velocity, magnetic field strength, field 
line curvature radii etc, which prevail in the radio emission region. For estimating the radio 
emission altitudes two kinds of methods have been proposed: (1) Purely geometric method, 
which assumes the pulse edge is emitted from last open field lines (e.g., Cordes 1978; Gil 
& Kijak 1993; Kijak & Gil 2003), (2) Relativistic phase shift method, which assumes the 
asymmetry in the conal components phase location, relative to core, is due to the aberration- 
retardation phase shift (e.g., Gangadhara & Gupta 2001, hereafter GGOl). By estimating the 
phase lag of polarization angle inflection point with respect to the centroid of the intensity 
pulse, Blaskiewicz et al. (1991) have estimated the emission heights. The results of purely 
geometric method are found to be in rough agreement with those of Blaskiewicz et al. (1991). 
However, compared to geometric method, the emission heights estimated from relativistic 
phase shift are found to be notably larger, particularly in the case of nearly aligned rotators 
(Gupta & Gangadhara 2003, hereafter GG03). Dyks, Rudak and Harding (2004, hereafter 
DRH04) by revising the relation for aberration phase shift given by GGOl, have re-estimated 
the emission heights. In the small angle approximation, they have found that the revision 
furnishes a method for estimating radio emission altitudes, which is free of polarization 
measurements and does not depend on a the magnetic axis inclination angle and (3 the sight 
line impact angle. By assuming the beamed radio waves are emitted in the direction of 
field line tangents, Gangadhara (2004, hereafter G04) has solved the viewing geometry in an 
inclined and rotating dipole magnetic field. 

In §2, we derive the angle between the corotation velocity of particles/plasma bunches 
and dipolar magnetic field. Using the magnetic colatitude and azimuth of the emission 
spot in an inclined and rotating dipole, we have derived the phase shift due to aberration, 
retardation and polar cap current in §3, and the emission radius in §4. In §5, we compare the 
shifts due to various process. As an application of our model, in §6, we have re-estimated 
the emission heights of six classical pulsars. Further, by considering a mean profile of a 
millisecond pulsar PSR J0437-4715, we have estimated the relativistic phase shift in conal 
components and their emission heights. 
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2. ANGLE BETWEEN PLASMA ROTATION VELOCITY AND DIPOLAR 

MAGNETIC FIELD 

Consider an inclined and rotating magnetic dipole (mt) with the rotation axis Cl, as 
shown in Figure 1. The angles a and (p' are the magnetic axis inclination angle and the 
rotation phase, respectively. Assume that the relativistic secondary plasma flows along the 
dipolar field lines, and emit the beamed radiation in the direction of field line tangents (6ot)- 
In a non-rotating case, to receive such an emission the sight line (n) must line up with &ot- 
Let Q be the emission point on a field line at which n = {sin(, 0, cos^) is parallel feot, where 
C = a + P and (3 is the sight line impact angle relative to rht- The angles ZQOR = 9 and 
ZQRT = are the magnetic colatitude and azimuth of Q relative to rht, respectively. While 
ZQOZ = 9' and ZXOS = 4>' are the colatitude and azimuth of Q relative to rotation axis, 
respectively. 

The position vector of Q is given by (sec cq. 2 in G04). For brevity we shall drop 
the suffix, and take r = r^t- The angle between r and Cl is given by 

cos 9' — Cl.r — cos a cos ^ — sin a sin 9 cos , (1) 

where the unit vector r = r/|r|, and the expressions for 9 and 4> as functions of 0', a and /3 
are given in G04. 

For the field line which lies in the meridional plane, defined by tl and rht, we have 
^min = + ^min- The magnetic colatitude 6'inin of Q can be obtained by setting 0' = in 
equation (9) of G04: 



/ cos(2/3) + V2cos/3v/17 + cos(2/3) 
' S 



1 / cosf2/3) + V2cos/3A/l7 + cosf2/3) - 1 
fmin = -arccosi g 



^ '^/3 + 0{/3f for P<^1. (2) 

This is the minimum value, which 9 takes at (j)' — 0. 

The rotation velocity of the plasma particle (bunch) at Q is given by 

Vrot = O, X r = Qrsm9' 6 , (3) 

where Q is the pulsar angular velocity, and the unit vector e represents the direction of 
rotation. Consider a Cartesian coordinate system - XYZ, with Z-axis parallel to the rotation 
axis and X-axis lies in the fiducial plane defined by h and fl. Let G be the angle between 
the field line tangent 6ot and e, then we have 

e = cos © e|| + sin © €± , (4) 
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where the unit vectors ey and €± are parallel and perpendicular to bot- Therefore, the angle 
© is given by 

cos = e.^ot = — = , (5) 

where 

Oi = sin q; sin (f), 02 = cot ^, 
cos^ a cos^ + a| sin^ a + 02 cos sin(2 a) + sin^ 



^3 - 1,2 

1 + ai 

We have plotted © as a function of rotation phase (j)' for different a in Figure 2. It shows 
© = 90° for the field lines which lie in the meridional plane {(f)' — 0) but for other field lines 
it is < 90° on leading side and > 90° in the traihng side. But for an ahgned rotator {a — 0°), 
it is 90° for all the field lines. 



3. PHASE SHIFT OF RADIATION EMITTED BY A PARTICLE (BUNCH) 

Since the pulsar spin rate is quite high, the rotation effects such as the aberration and 
retardation play an important role in the morphology of pulse profiles. For an observer in 
an inertial frame, the radiation beam gets phase shifted due to the corotation of plasma 
particles and the difference in emission radii. 

Since the radiation by a relativistic particle is beamed in the direction of velocity, to 
receive it sight line must align with the particle velocity within the angle I/7, where 7 is the 
Lorentz factor. The particle velocity is given by 

V = Kcfeot + Vrot, (6) 

where c is the speed of light. By substituting for Vrot from equation (3) into equation (6) we 
obtain 

V = {kc + fir sin 9' cos ©) 6ot + sin 9' sin © e± . (7) 
By assuming |v| ~ c, from equation (7) we obtain the parameter 



/ Or \ ^ Or 
K= ^jl - l—j sin^ ^' sin^ © - — sin ^' cos © . (8) 

In Figure 3, we have plotted k as a function of r for different a. It shows k ~ 1 for r/r^c ^ 1, 
but at large r it decreases from unity due to increase in rotation velocity, where r^c is the 
hght cyhnder radius. Machabeh & Rogava (1994), by considering the motion of a bead inside 
a rotating linear tube, have deduced a similar behavior in the velocity components of bead. 
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Using equation (3) we can solve equation (4) for e±, and obtain 



= ■ n, ■ - cot © &0t ■ (9) 

sm d' sm c) 



Let %l) be the angle between the rotation axis and v, then we have 



cosijj — Ct.v — cos( I sin^ ^' sin^ © — — sin ^' cos © 

= K cos( , (10) 
where v — v/|v|. For r/r^c <^ 1, it reduces to ~ C- 

3.1. Aberration Angle 

If 7] is the aberration angle, then we have 

- ^ + i7r sin 6*' cos © , , 

COS77 = bofV^ , (llj 

sin 77 = e± • ■?) = -j — j- sin ^' sin©. (12) 

Therefore, from equations (11) and (12), we obtain 

0,r sin ^' sin© 

tan77 = , (13) 

c ^1 - (l]r/c)2 sin^ ^'sin^© 

Hence the radiation beam, which is centered on the direction of v, gets tilted (aberrated) 
with respect to 6ot due to rotation. 

For flr/c <^ 1, it can be approximated as 

Qr 

tan T] a; — sin 9' sin © . (14) 
c 



3.2. Aberration Phase Shift 

Consider Figure 4 in which ZAD, ZBX, ZCY and DXY are the great circles centered on 
the neutron star. The small circle ABC is parallel to the equatorial great circle DXY. The 
unit vector hot represents a field line tangent, which makes the angle ( with respect to the 
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rotation axis ZO. The velocity unit vector v is inclined by the angles rj and t/j with respect 
to 6ot s-nd ZO, respectively. We resolve the vectors 6ot ^-nd v into the components parallel 
and perpendicular to the rotation axis: 

6ot = sin C 6ot± + cosC fl , (15) 

■u = sin-^ v± + cosip Cl , (16) 

where the unit vectors 6ot± and lie in the plane of small circle ABC. Next, by solving for 
6ot± and v±, we obtain 

&ot± = -^(^ot-cosC^^) , (17) 
smC 

v± — — ^ — (-0 — cosip Cl ). (18) 

sin ip 

By taking scalar product with 6ot± on both sides of equation (18), we obtain 

cos((50^be) = ^-L • ^ot± = ■ bot± - COSl/j Cl ■ bot± ). (19) 

sm-^ 

Since Cl and 6ot± are orthogonal, we have 

cos(50lbe) = -i-(t)-6ot±). (20) 

smyj 

Using bot± from equation (17) we obtain 

X 1 {v - hot- cosC V - Cl) 

cos 50^be) = — — ■ 21) 

sm ip sm 

By substituting for v ■ 6ot and v • Cl from equations (11) and (10), we obtain 

_ ,, , 1 (cos 77 — cos C cosih) 

COS ^^^be = ^-T —. 22 

^ ^"^^ smV' smC 

Substituting for rj again from equation (11), we obtain 

1 f fir \ 

cos(50' be) = { K-\ sin ^' cos © — cos C cos ) . (23) 

sm^sm-^ \ c ) 

By substituting for k from equation (10), we obtain 

1 / cos flv \ 

cos(50' Up) = H sin Q' cos © — cos C cos ) ■ (24) 

smCsm'^ V^osC c ) 
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It can be further reduced to 

/r,/ \ > , ilrsin^'cos© 

cos(d0abe) = tan C cot + — -. (25) 

c sm Q sm ip 

Next, by taking a scalar product with tl x 6ot± on both sides of equation (18), we obtain 

sin((5(/);j ^ (Q X 6ot±) • ^± = ■ A (lxhot)-v . (26) 

sm Q sm 

By substituting for v from equation (7), we obtain 

. , Or sin 6*' sin 9, - - , ^ , „, 

sm<50lbeH — r-— O^xfcot -ex ■ 27 

c sm (, sm ip 

Using e±_ from equation (9), it can be further reduced to 

Or A 



where 



A = (Q X 6ot) • (^^ X ^)/sinC 

= cos 0'(cos ^ sin a + cos a cos (j) sin 6') — sin (f)' sin sin ^. 

Thus, we obtain from equations (25) and (28): 

, , , (f2r/c) sin 

tan( 0^be) = sin ( tan C cos ^ + (fir /c) sin ^' cos 9 ' ^ ^ 

The aberration phase shift ^^^bc plotted as a function of /9 in Figure 5 for different a in 
the two cases: (1) 0' = 0° and r /r-^^, = 0.1, and (2) 0' = 60° and r/r^^ = 0.1. The Figure 5a 
shows for a ~ 90°, 50abc ~ ''"/''"lc' which is nearly independent of f3. But, for other values 
of a, it does depend on P : larger for /? < and smaller for P > 0. Since 6'' decreases with 
10' I, and 9 < 90° on leading side and > 90° on traihng side, 50abe smaller in Figure 5b 
compared to its corresponding values in Figure 5a. Further, we note that 50abe negative 
gradient with respect to |0'| for both the signs of p. In Figure 6a, 50abe plotted as a 
function of 0' for different a and fixed /3 = 4°. It is highest at large a and small |0'|. 

For r/rLc ^ 1, we can series expand 50abc ^^^^ obtain 

Scl>',^, = hr^ + hrl + 0{r^f , (30) 

where rn = r/rLc, 

smC 

b2 — — cot^ C COS 9 sin 9' bi . 
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3.3. Retardation phase shift 

Let 9c, 0c and fc be the magnetic colatitude, azimuth and the position vector of the 
emission spot at the emission time. In the expressions for 6 and (see eqs. 9 and 11 in G04) 
we replace (p' by (p' + ^^^be obtain 6'e and 0e- Then using the values of 9^ and 0e we find 
the unit vector fe (see eq. 2 in G04). For brevity, we shall drop the suffix on f. Note that 
0' > on the trailing side (Fig. 1). Consider the emission radii ri and r2 such that ri < r2. 
The time taken by the signal emitted at the radius ri is given by 

ti — -(d — ri • n), (31) 
c 

where d is the distance to the pulsar, and n — (sin^, 0, cos^) is the unit vector pointing 
toward the observer. For another radius r2, the propagation time is given by 

t2^-{d-r2-n), (32) 
c 

The radiation emitted at the lower radius ri takes more time to reach the observer than the 
one emitted at r2. The time delay between the two signals is given by 

St = ti — t2 = -(r2 ■ n — ri ■ fi). (33) 
c 

By considering the neutron star center (ri = 0) as the reference and r2 = r, we obtain 

5t^^{f-h). (34) 

Let (7 be the angle between f and n, then we have 

cosa = f-n = cosC(cosQ;cos^e — cos^esinasin^e) — sinC[sin0'sin0eSin6'e — 

cos 0'(cos 9e sin a + cos a cos 0e sin ^e)]- (35) 

The time delay 6t of components emitted at lower heights, shifts them to later phases 
of the profile by (e.g., Phillips 1992) 

^^'ret — = COSCT . (36) 

In Figure 6b, 50ret is plotted as a function of 0' for different a and fixed /3 = 4°. At small 
a, it is nearly constant and larger. But at larger a it falls with respect to |0'|, as f inclines 
more from n. 

For Tn 1, we can find the series expansion: 

5(f)[,, = ci rn + C2 rl + 0{r^^ , (37) 
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where 

ci = cos(r — 6) , 

3 sin a sin^ sin(2 9) sm(T ~ ^) ^ 
2 sin C sin 0'V8 + cos^ V 

and r is the half-opening angle of the emission beam (see eq. 7 in G04) . 



3.4. Relativistic Phase Shift 

Since the retardation and aberration phase shifts are additive, they can collectively 
introduce an asymmetry into the pulse profile (e.g., GGOl). Therefore, the relativistic phase 
shift is given by 

= rnCOSCT + arctan ( ^ — — rnSinCA ^\ 

\sm(; tanC cosi/y + rnSm 0' cosB / 

In Figure 7, we have plotted ^^^pg as a function of 0' for different a in the two cases of 
P — ±4°. It shows 50rps reaches maximum at 0' ~ 0, and falls at large 

In the limit of Tn <C 1, we can series expand 50rps obtain 

50rps = A*i + A«2 rl + 0{rnf , (39) 

where 

/xi = cos(r -9) + - 



sinC ' 

/3sina sinCsin^0sin(2 6') sinfr — 6*) 9^ ^ ^/A ^ 
112 = ^ — ^ - cos^Ccos©sm6'' ' 



2 sin0V8 + cos2r / sin^C 



3.5. Phase Shift Due to Polar Cap Current 

Goldreich and Julian (1969) have elucidated that the charged particles relativistically 
stream out along the magnetic field lines of neutron star with aligned magnetic moment and 
rotation axis. Hibschman and Arons (2001) have shown that the field-aligned currents can 
produce perturbation magnetic field Bi over the unperturbed dipole field Bq, and it can 
cause a shift in the polarization angle sweep. Here we intend to estimate the phase shift in 
the intensity profile due to the perturbation field Bi. We assume that the observed radiation 
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is emitted in the direction of tangent to the field B = Bq + Bi. Using the equations D5 and 
D6 given by Hibschman and Arons (2001), we find the Cartesian components of perturbation 
field: 

cos a sin sin 



Bi = 



H cos q; sin ^ cos 
^ , 



(40) 



where is the magnetic moment. The Cartesian components of unperturbed dipole is given 

by 

3 11 



Bn 



2r 



3 / 

■sin(2^)cos0, -4sin(2^)sin0, 4(3cos^^-l) 

2 TfA J.6 



The magnetic field, which is tilted and rotated, is given by 

Bt = A B, 

where A is the transformation matrix given in G04. 

The component of Bt perpendicular to the rotation axis is given by 

Btx = Bt - (Bt.Q)Q. 
Similarly, we find the perpendicular components of Bgt = A • Bq : 

Bot± — Bot — (Bot-f^)f^- 
If 50pp is the phase shift in Bt due to the polar cap current then we have 



(41) 



(42) 



cos ((50' 



(43) 
(44) 

(45) 



where ht± = Btx/|Bt^|, the unit vector 6ot± = Bot±/|Bot±| is parallel to the unit vector x 
along the X-axis, and -Bt±,x is the X-component of Btj,. If y is the unit vector along Y-axis, 
then we have 



B 



t±\ 



(46) 



Therefore, we have 



tan((5(/)[.J 



t-L,y 



where 



-Bt±,x 

C?2 + 4 '^n ' 



di — cos0'sin(2Q;) — 2cos^Q;tanC , 
d2 — 3 cos ^ tan C , 
ds — — sin(2Q;) sin0' . 



(47) 
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We have plotted 54)'^^ as a function of 0' in Figure 8 at a = 10°, r-^ — 0.01 and 0.1. It 
decreases with the increasing \<p'\ and it is mostly negative, except in the case of /3 < 0, 
where it is positive over a small range of 0' near rht) plane. So, tries to reduce the 
relativistic phase shift (^^rpg, except over a small range of <p' where it enhances the shift in 
the case of /3 < 0. 

In the limit of Tn <S 1, we can series expand equation (47) and obtain 

'^<^;c = |rn-^r^ + 0(r„)l (48) 



4. EMISSION RADIUS FROM PHASE SHIFT 

We can find the net phase shift due to aberration, retardation and polar cap current by 
adding equations (38) and (47): 

/ Tn sin ^ A \ 

= Tn cos (T + arctan - — — ■ ■. — — + 

\sm Q tan Q cos ip -\-r^ sm U' cos fc) / 



arctan 



d2 + Tr, 



(49) 



In Figure 9 we have plotted 50' as a function of 0' in the four cases of (0.01, 0.1, 0.2 and 
0.3). It shows 50' reaches maximum near the {Q, rht) plane and falls with respect to |0'|. 
Note that the magnitude of gradient of 6(f)' with respect to |0'| is higher than that of 50rps- 
In the case of /3 < 0, we find 50' becomes negative at large |0'| as the magnitude of 50;^, 
exceeds the magnitude of 5(f)[.pg. At higher rn, we note that 50' is slightly asymmetric about 
0' = O. 

For Tn <^ 1, we obtain 

50' = z/lr„ + z/2r^ + O(r„)^ (50) 
where I'l = + {di/d2) and 1^2 = 1^2 — {di d^/dl). 

For 50' <^ 1, we can solve equation (50) for the emission radius and obtain 

r^'J^ScI)' -"^^Scl)'^ + 0{S<P'f. (51) 



1 
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5. DISCUSSION 

There are many processes such as aberration, retardation, polar cap currents etc, which 
can collectively introduce a phase shift in the pulse components. To compare the magnitude 
of shifts due to each one of these processes, we estimate the shifts in the order of r^. 



5.1. Relativistic phase shift 50rps in the order of r^ 

To estimate the order of magnitude of in the order of rn, we set f^^rpg = rn^% and 
obtain 

U = (52) 

In(rn) 

We have plotted ^rps as a function of 0' in Figure 10 for rn = 0.01 and 0.1. It is less than 
unity in both the cases. 



5.2. Magnitude of refinement in in the order of r-^ 

In the limit of small angle {9' ~ C,), i.e., in the emission region close to the meridional 
(n, rhi) plane and for rn <S 1, it can be shown that our expression for relativistic phase 
shift (eq. 38) reduces to r ^ (^lc/2)50o7^^ given in DRH04, where S(f)']^jiff is the relativistic 
phase shift. We can estimate (^0difr,i — '^^rps ~ ^'P'drh^ ^■^■■> difference in the phase shifts 
predicted by the two formulas. By substituting for (50^pg from equation (38) and i^^o 
the expression for (50diff,n series expanding it in powers of rn, we obtain 

H'm.i = (/^i - 2)r„ + ^X2rl + O(r^). (53) 

In Figure 11, we have plotted the factor (//i — 2) which appears in the leading term of above 
equation as a function of phase 0' for the two cases: (i) a — 10° and (3 — —4°, 0°, 4° and 
(ii) a = 90° and /3 = 0°, ± 20°, ± 40°. It shows the factor (/xi - 2) varies from zero with 
respect to 0' as well as j3. At 0' = and /9 = 0, i.e., along the magnetic axis in the meridional 
(O, rhi) plane we find — 2, and hence (^0diff,i becomes third order in rn. 

To estimate the magnitude of the leading term in equation (53) in the order of rn, we 
set {jii — 2)rn = r^'^'^ , and obtain the index 

hi(|//i-2|) 

4difr,i = 1 + — T-j—. — , 54 
m(rn) 

where we consider the absolute values of {jii — 2) as jii can have values > as well as < 2. 
In Figure 12 we have plotted ^difr,i as a function of 0' in the cases of different a and (3. The 
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panels (a) and (c) for rn = 0.01 shows ^difr,i < 2 at almost all the longitude except at the 
spikes where fj^i — 2 crosses the zero. The panels (b) and (d) for larger = 0.1, shows ^diff,i 
is larger than 2 for \/3\ < 1 near the {Cl, rht) plane, but at large \(f)'\ it lies between 1 and 2. 



5.3. Phase shift due to Polar cap current Scj)'^^ in the order r^ 



To estimate the magnitude of Scf)^^ in the order rn, we define Scp^^ = {rnY^", and obtain 

e.^!^^. (55) 

ln(r„) 

Since 50p^ is mostly negative except in the case of (3 < where it is positive near the 
meridional plane, we use the absolute values of S(f)p^ in computing ,^pc. We have plotted 
as a function of 0' in Figure 13. It shows is of the order of 3/2 near 0' ~ and unity at 
large \(f)'\ except at the spikes. 



5.4. Net phase shift 5(1)' in the order rn 
To estimate S4>' in the order rn, we define — r^ and obtain 

(56) 

ln(r„) 

In Figure 14 we have plotted ^ as a function of 0'. It is of the order of unity in the region 
close to 0' ~ and greater than unity at large 

Again, we can estimate difference in phase shifts 50diff,2 — ^4>' ~ ^4>'drh predicted by the 
two formulas. If Scp'^^Q ^ — f^^^^^ then we have 

In(rn) 

In Figure 15 we have plotted ^diff,2 as a function of 0'. It shows, except in the spiky regions, 
Cdifr,2 is ~ 3/2 near 0' ~ and < 1 at large 



5.5. Magnetic field sweep hack 



Due to the rotational distortions such as the magnetic field sweep back of the vacuum 
dipole magnetic field hues, the relativistic phase shift is likely to be reduced. The magnetic 
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field sweep back was first considered in detail by Shitov (1983). Recently, Dyks and Hard- 
ing (2004) investigated the rotational distortions of pulsar magnetic field by assuming the 
approximation of vacuum magnetosphere. We used their expressions (eqs. (12) and (13) in 
Dyks and Harding, 2004), to estimate the magnetic field sweep back: 

^2 



1/ 



M^-t 2 



XZ I ^ \ 

3— cos a + 3— — 1 sm a 



K- (58) 



Using X = r sin 9' cos (f)' and z = rcos9', we have plotted 50mfsb ^ function of in 
Figure 16a using 0' = 50°, (3 = 0°, a = 10° and 90°. It shows S(f)[^fg^ increases with Tn, 
and higher in the case of orthogonal rotators. However, it is smaller than the aberration, 
retardation and polar cap current phase shifts for rn < 0.2. By defining 50mfsb = r^^f^b^ 
obtained 

U.. = if^. (59) 
In(rn) 

In Figure 16b, we have plotted .^mfsb as a function of r^. It shows the index ^mfsb > 3 in the 
case of a = 10°, while in the case of o; = 90° it lies in the range of 2 < ^mfsb < 3. 

In addition to the various processes, which have been considered, the corotation of 
Goldreich- Julian charge density {i]qj) can also produce the phase shift. The corotating 
charges induces magnetic field B^ot given by 

47r 

V X Brot = — 77GJ X r. (60) 
c 

By defining Brot = Brot,i + Brot,2, and the expression for f2 x r from equation (3), we can 
resolve equation (60) into two component equations: 

4:71 

V X Brot,i = —Vgj ^ r sin 9' sin 9 e± , (61) 

4:71 

VxBrot2 = — rjG^^lr sin 9' cos Qe\\ . (62) 

c 

The equation (61) implies 6ot ■ V x Bfot,! = 0. Since Bi-ot,i lies almost parallel to 6ot in the 
emission region, we do not expect any significant phase shift due to Brot,i- But Brot, 2 can 
introduce a phase shift. Since both flr/c and rrjcj/Bo are first order in r^, -Brot,2/-Bo becomes 
second order in rn. In addition © ~ 90°, therefore, we neglect the phase shift due to Brot,2- 

Finally, we may summarize that among the various phase shifts considered the relativis- 
tic phase shift due to aberration-retardation is the dominant, as indicated by the Figure 10. 
In the small angle approximation, i.e., when the range of 4>' and /3 are small, our equation (38) 
for the relativistic phase shift reduces to the expression given by DRH04. The neglected ef- 
fects such as the magnetic field fines sweep back due to the reaction force exerted by the 
magnetic dipole radiation and the toroidal current due to the corotation of magnetosphere 
are of higher order than the proposed refinement. 
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6. APPLICATION 

The conal components are believed to arise from the nested hollow cones of emission 
(Rankin 1983a,b, 1990, 1993), which along with the central core emission, make up the pulsar 
emission beam. The observed phase shift of conal pair with respect to core is given by 

S<Po^. = 5</'cone - ^0cr , ^'/'cone > <^0cr > 0, (63) 

where Scj)'^^^^ is the phase shift of cone center and is the phase shift of core peak with 
respect to the meridional plane. In the relativistic phase shift model, we assume that the 
core is emitted from the lower altitudes and hence 50^^ ~ as the aberration and retardation 
effects are minimal at those heights. 

Let be the peak locations of conal components on leading and trailing sides 

of a pulse profile, respectively. Then, using the following equations (GGOl), we estimate the 
phase shift of cone center with respect to core and the phase location of component peaks 
in the absence of phase shift, i.e., in the corotating frame: 

SK^s = -^(0T + 0l) , 0' = ^(0T - 0l) ■ (64) 



6.1. Classical pulsars 

The emission height rem of PSR B0329+54 has been estimated by GGOl and six more 
pulsars by GG03. Later, DRH04 have revised the aberration formula of GGOl and re- 
estimated the emission heights of all those seven pulsars. Using the revised phase shift 
given by equation (49) and S4>' ~ 50obs' *® have computed the emission heights of all those 
pulsars considered in GGOl and GG03, except PSR B2045-16. Since the emission heights 
of PSR B2045-16 reported in GG03 and DRH04 are notoriously large, we have dropped it 
from this study. The drifting phenomenon of PSR B2045-16 (e.g., Oster & Sieber 1977) 
has probably complicated the identification of the component peak locations. In Table 1, 
we have given the cone number in column 4 and the revised emission heights in columns 5. 
Note that these emission heights are measured upwards from the core emission region. In 
the present calculations, we have assumed core emission height is zero. However, if there is 
any finite emission height for the core then the emission heights of all the cones will increase 
correspondingly. To compare the emission heights predicted by our formula (eq. 49) with 
those by equation (7) in DRH04, we define the percentage of refinement A as 

A = (^em-/^L) ^QQ^ ^gg^ 
"'em 
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Table 1. Radio emission altitudes /lem and refinement A. The numbers are based on the 
cone shifts measured by GGOl and GG03, and a. and /3 values from Rankin (1993). 



Pulsar 


P[s] 


V [MHz] 


Cone"- 


feem [km] 


^■em [% of TLC] 


A[%] 


sj Slof 


B0329+54 


0.7145 


325 


1 


154±80 


0.45 


3.12 


0.58±0.15 






325 


2 


339±59 


0.99 


3.39 


0.56±0.05 






325 


3 


619±84 


1.82 


3.84 


0.56±0.04 






325 


4 


921±250 


2.70 


4.73 


0.64±0.09 






606 


1 


123±77 


0.36 


3.11 


0.64±0.20 






606 


2 


292±55 


0.86 


3.34 


0.57±0.06 






606 


3 


479±140 


1.40 


3.71 


0.60±0.09 






606 


4 


670±180 


1.96 


4.46 


0.69±0.10 


B0450-18 


0.5489 


318 


1 


248±18 


0.95 


7.10 


0.63±0.02 


B1237+25 


1.3824 


318 


1 


161±40 


0.24 


0.01 


0.35±0.04 






318 


2 


415±29 


0.63 


0.06 


0.43±0.02 






318 


3 


536±23 


0.81 


0.15 


0.62±0.01 


B1821+05 


0.7529 


318 


1 


235±100 


0.65 


2.43 


0.48±0.10 






318 


2 


332±94 


0.92 


2.76 


0.64±0.09 






318 


3 


456±87 


1.27 


3.15 


0.71±0.07 


B1857-26 


0.6122 


318 


1 


170±53 


0.58 


4.31 


0.70±0.11 






318 


2 


372±30 


1.27 


5.50 


0.78±0.03 






318 


3 


536±82 


1.83 


6.72 


0.84it0.06 


B211H-46 


1.0147 


333 


1 


878±210 


1.81 


8.55 


0.27±0.03 






333 


2 


1420±54 


2.94 


13.60 


0.38±0.01 



Cone numbering is the same as in GGOl and GG03 (i.e. from the innermost cone outwards). 
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where hem is the emission height estimated from equation (49), and h'^^ from the equation (7) 
in DRH04. In column 7, we have given the values of A. We note that the refinement increases 
from inner cone to outer cones for any given pulsar in the table. It is least (< 1%) in the 
case of PSR 1237+25 but greater than 2% for all other pulsars. It is maximum (13%) in the 
case of PSR B2111+46. In column 8, we have given the colatitude s of foot field lines (eq. 15 
in GGOl) relative to magnetic axis on the polar cap. It is normalized with the colatitude 
siof of last open field line. 

6.2. Millisecond pulsar: PSR J0437-4715 

We consider the nearest bright millisecond pulsar PSR J0457-4715, which has a period 
of 5.75 ms. It was discovered in Parkcs southern survey (Johnston et al. 1993), and is in 
a binary system with white dwarf companion. Manchester & Johnston (1995) have pre- 
sented the mean pulse polarization properties of PSR J0437-4715 at 1440 MHz. There is 
significant linear and circular polarization across the pulse with rapid changes near the pulse 
center. The position angle has a complex swing across the pulse, which is not well fitted with 
the rotating vector model, probably, due to the presence of orthogonal polarization modes. 
The mean intensity pulse has a strong peak near the center, where the circular polarization 
shows a clear sense reversal and the polarization angle has a rapid sweep. These two features 
strongly indicate that it is a core. The profile shows more than 8 identifiable components. 



Consider the average pulse profile given in Figure 17 for PSR J0437-4715 at 1440 MHz. 
To identify pulse components and to estimate their peak locations wc followed the procedure 
of Gaussian fitting to pulse components. Two approaches have been developed and used by 
different authors. Unhke Kramer et al. (1994), who fitted the sum of Gaussians to the total 
pulse profile, we separately fitted a single Gaussian to each of the pulse components. We used 
the package Statistics'NonhnearFit' in Mathematica (version 4.1) for fitting Gaussians to 



Table 2. Parameters related to radio emission from PSR J0437-4715 at 1440 MHz 



Cone 




<I>'t 




r 






A 






n 


n 


n 


n 


(Km) 


(%r-Lc) 


(%) 




3 


-72.00±0.15 


47.44±0.08 


12.30±0.08 


18.00it0.02 


35.90±0.25 


13.10 


18.00 


0.57±0.00 


4 


-92.51±0.36 


65.64±0.29 


13.40±0.23 


22.90±0.05 


49.90±0.87 


18.20 


35.50 


0.61±0.01 


5 


-120.00±1.01 


99.04±4.15 


10.50±2.10 


29.30±0.40 


63.40±9.60 


23.10 


60.40 


0.69±0.05 
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pulse components. It can fit the data to the model with the named variables and parameters, 
return the model evaluated at the parameter estimates achieving the leastsquares fit. The 
steps are as follows: (i) Fitted a Gaussian to the core component (VI), (ii) Subtracted the 
core fitted Gaussian from the data (raw), (iii) The residual data was then fitted for the next 
strongest peak, i.e., the component IV, (iv) Added the two Gaussians and subtracted from 
the raw data, (v) Next, fitted a Gaussian to the strongest peak (IX) in the residual data. 
The procedure was repeated for other peaks till the residual data has no prominent peak 
above the off pulse noise level. 

By this procedure we have been able to identify 11 of its emission components: I, II, 
III, IV, V, VI, VII, VIII, IX, X and XI, as indicated by the 11 Gaussians in Figure 17. The 
distribution of conal components about the core (VI) reflect the core-cone structure of the 
emission beam. We propose that they can be paired into 5 nested cones with core at the 
center. In column 1 of Table 2 we have given the cone numbers. Since the assumption that 
^'core -C Tcone is uot fulfilled for the inner cones (1 and 2), we have not included their results. 
The peak locations of conal components on leading and trailing sides are given in columns 2 
and 3, respectively. In column 4 we have given the values of 50obs- 

Manchester & Johnston (1995) have fitted the rotating vector model (RVM) of Rad- 
hakrishnan & Cooke (1969) to polarization angle data of PSR J0437-4715 and estimated 
the polarization parameters a = 145° and ( = 140°. However, for these values of a and ( 
the colatitude s/siof (eq. 15 in GGOl) exceeds 1 for all the cones. The parameter s gives 
the foot location of field lines, which are associated with the conal emissions, relative to 
magnetic axis on the polar cap. It is normalized with the colatitude siof of last open field 
line. However, Gil and Krawczyk (1997) have deduced a — 20° and /3 = —4° by fitting 
the average pulse profile of PSR J0437-4715 at 1.4 GHz rather than by formally fitting the 
position angle curve. Further, they find the relativistic RVM by Blaskiewicz et al. (1991) 
calculated with these a and /5 seems to fit the observed position angle quite well than the 
non-relativistic RVM. We have adopted these values of a and /3 in our model as they also 
confirm s/siof < 1 (see below). 

The half-opening angle F of the emission beam is given in column 5, and it is about 30° 
for the last cone. Using the revised phase shift given by equation (49) and S(p' ~ 50^^^, we 
have computed the emission heights given in column 6, and their percentage values in tlc 
in column 7. It shows that the third cone is emitted at the height of about 13 Km and all 
other successive cones are emitted from the higher altitudes. Note that these conal emission 
heights are measured upwards from the core emission region. In the present calculations, 
we have assumed that the core emission height is approximately zero. However, if there is 
any finite emission height for the core then the emission heights of all the cones will increase 
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correspondingly. Based on the purely geometric method Gil and Krawczyk (1997) have 
estimated the emission height of 44 Km. Of course, their method is a symmetrical model 
which assumes the pulse edge is emitted from the last open field lines. 

The refinement A values are given in column 8. It is about 18% at the third cone and 
60% at the outer most. In column 9, we have given the colatitude s/siof of foot field lines, 
which are associated with the conal emissions. It shows due to the relativistic beaming and 
geometric restrictions, observer tends to receive the emissions from the open field lines which 
are located in the colatitude range of 0.6 to 0.7 on the polar cap. 

7. CONCLUSION 

We have derived a relation for the aberration phase shift, which is valid for the full range 
of pulse phase. Though in the small angle approximation we can show that the aberration 
phase shift becomes independent of parameters a and it docs depend a and /? in the 
case emissions from large rotation phases or altitudes. We have given the revised phase 
shift relation by taking into account of aberration, retardation and polar cap current. We 
find among the various phase shifts considered the relativistic phase shift due to aberration- 
retardation is the dominant. 

The emission heights of six classical pulsars have been recomputed, and analyzed the 
profile of a millisecond pulsar PSR J0437-4715. In the profile of PSR J0437-4715 we have 
identified 11 of its emission components. We propose that they form a emission beam 
consist of 5 nested cones with core component at the center. The emission height increases 
successively from cone 3 to cone 5. 

In the limit of small angle [fi ~ 0) and low altitude rn <^ 1) approximation, our 
expression for relativistic phase shift (eq. 49) reduces to equation (7) given by DRH04. In 
the case of pulsars with small /3 and narrow profiles, we can use the approximate expression 
(eq. 7) given in DRH04 to estimate the emission heights while for the pulsars with wide 
profiles or large /3 we have to use the revised phase shift given by our equation (49). 

I thank R. M. C. Thomas, Y. Gupta for discussions, and J. A. Gil and J. Murthy for 
comments. I am thankful to M. N. Manchester and S. Johnston for providing EPN data, 
and the anonymous referee for useful comments. 
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Fig. 1. — Viewing geometry of emission beam. The heavy hne elhpse represents the cone of 
emission centered on the magnetic axis rhf The arcs ZQX, ZRS, ZY and XSY represent the 
great circles centered at O (star center). The magnetic colatitude and the phase angle (f)' 
of the emission spot are measured from the meridional (Jl, TTit) plane. They have the signs 
in such a way that 0' is positive while is negative on the trailing side, and the vice versa 
on leading side. 
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Fig. 2. — The angle © vs rotation phase 0' is plotted for different values of a and a fixed 
/3 = 2°. 
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Fig. 4. — Celestial sphere describing the aberration phase shift of pulsar radio emission, 
where rj is the aberration angle, and 50abe ^he corresponding phase shift. 
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Fig. 5. — Aberration phase shift ^^^bc vs the sight hne impact angle j3 for different a : panel 
(a) for Tn = 0.1 and 0' = 0°, and (b) for = 0.1 and 0' = 60°. 




Fig. 6. — Aberration and retardation phase shifts vs the phase (j)' . Chosen — 0.1, /3 = 4° 
and different a for both panels. 
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Fig. 8. — The phase shift 5(f)'^^ due to polar cap current vs phase 0'. Chosen a — 10° and 
rn = 0.01 for panel (a), and rn = 0.1 for (b). 
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Fig. 9. — The net phase shift due to aberration, retardation and polar cap current 54>' vs 
phase (f)'. Chosen a — 10° and Tn = 0.01, 0.1, 0.2 and 0.3 for panels (a), (b), (c) and (d), 
respectively. 




Fig. 10. — The index ^rps vs phase 0'. Chosen a = 10° and r-a — 0.01 for panel (a), and 
a = 10° and Tn = 0.1 for (b). 
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Fig. 11. — The factor (//i — 2) vs the phase 0' at a fixed a and different /3. The panel (a) is 
plotted using a = 10° and /3 = -4°, 0°, 4°, and (b) with a = 90° and f3 = 0°, ±20°, ±40°. 




Fig. 12. — The index ^diff,i vs phase (f)' at a fixed a and different p. Chosen for panel (a) 
rn = 0.01, a = 10°, (b) = 0.1, a = 10°, (c) = 0.01, a = 90°, and (d) Tn = 0.1, a = 90°. 
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Fig. 14. — The index ^ vs phase 0'. Chosen a — 10° and rn = 0.01 for panel (a), and a — 10° 
and Tn = 0.1 for (b). 
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Fig. 16. — The magnetic field sweep back 50mfsb and the index ^mfsb are plotted as functions 
of r/rLc (fn) in panels (a) and (b), respectively. Chosen 0' = 50°, P — 0°, a — 10° and 90°. 



-31 - 




Fig. 17.— Average pulse of PSR J0437-4715 at 1440 MHz and the model Gaussians (dotted 
line curves) fitted to emission components. 



